NASA— TM— 111579 


Plate mode velocities in graphite/epoxy plates 
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Measurements of* the velocities of the extensional and flexural plate modes were made along three 
directions of propagation in four graphite/epoxy composite plates. The acoustic signals were 
generated by simulated acoustic emission events (pencil lead breaks or Hsu-Neilsen sources) and 
detected by broadband ultrasonic transducers. The first arrival of the extensional plate mode, which 
is nondispersive at low frequencies, was measured at a number of different distances from the source 
along the propagation direction of interest. The velocity was determined by plotting the distance 
versus arrival time and computing its slope. Because of the large dispersion of the flexural mode, a 
Fourier phase velocity technique was used to characterize this mode. The velocity was measured up 
to a frequency of 160 kHz. Theoretical predictions of the velocities of these modes were also made 
and compared with experimental observations. Classical plate theory yielded good agreement with 
the measured extensional velocities. For predictions of the dispersion of the flexural mode, Mindiin 
plate theory, which includes the effects of shear deformation and rotatory inertia was shown to give 
better agreement with the experimental measurements. 

PACS numbers: 43.40.Dx 


INTRODUCTION 

Acoustic signals propagate in thin plates as the exten- 
sional and flexural plate modes when the wavelength of the 
acoustic signal is much larger than the plate thickness. This 
fact is of importance in a number of nondestructive evalua- 
tion (NDF) techniques that use acoustic waves such as the 
acoustic emission < AH) technique. The propagation of AH 
signals as plate modes has been demonstrated by a number 
of researchers. It was shown using simulated AE sources 
(pencil lead breaks) by Gorman 1 on thin aluminum and gr/ep 
composite plates and by Gorman and Prosser on thin alumi- 
num plates. A typical signal from a simulated AF source 
(pencil lead break) is shown in Fig. I. The extensional and 
flexural mode components of this waveform arc identified in 
this figure. AF signals from transverse matrix cracking 
sources in gr/ep composite plates were also shown to propa- 
gate as plate modes by Gorman and Ziola. 3 Smith 4 showed 
that crack growth events in thin aluminum plates under spec- 
trum fatigue loading produced signals that propagated as 
plate modes. 

The propagation of acoustic signals as plate modes in 
thin plates also occurs in other NDE techniques such as con- 
ventional ultrasonic testing, acousto-ultrasonics, and laser 
generated ultrasonics. Duke et r//. 5 suggested that the 
acousto-ultrasonic technique generated plate modes, made 
tone burst measurements of the velocities of the plate modes, 
and compared them with predictions based on classical plate 
theory (CPT). Schumacher et al. h demonstrated the existence 
of plate modes in laser generated ultrasonic signals in thin 
plates which were detected by noncontacting laser interfer- 
ometers. They also made measurements of the dispersion of 
these two modes in steel plates and compared them with 
theory and finite element modeling. 


For all of these NDE techniques, the propagation char- 
acteristics of these plate modes are of great importance. In 
acoustic emission testing, Gorman 1 discussed how erroneous 
source location could be obtained using conventional first 
threshold crossing or peak arrival techniques because of the 
presence of plate modes which propagate with different and 
dispersive velocities. Such source location errors were sub- 
stantiated by Ziola and Gorman and an alternative method 
for source location based on cross-correlation of the flexural 
mode waves was demonstrated. In the other mentioned NDE 
techniques, measurements of the velocities of the plate 
modes are useful for determining material properties and lo- 
cating flaws. Mai et af* Veidt and Sayir,'* and Dean 10 have 
discussed how measurements of flexural mode dispersion 
might be useful in determining the elastic constants of com- 
posite plates. 

In this research, measurements were made of the exten- 
sional and flexural plate mode velocities on four different 
composite laminates in three directions of propagation. The 
ply layups for these laminates were [0 U ,J, [0 4 ,90 4 | s , 
[0,W)| 4s , and ( 0,45, -45, ( J0] :> . Pencil lead breaks (Hsu- 
Neilsen sources) were used to generate the acoustic signals 
which were detected by wide band ultrasonic sensors. The 
first arrival time of the extensional mode as a function of 
distance from the source was used to compute its velocity. A 
Fourier phase technique was used to determine the flexural 
mode dispersion up to a frequency of 160 kHz. The measure- 
ments were made along the two principal axes (0 and 00 deg) 
in the plane of the plate and at an angle of 45 deg. 

The ability of existing plate theories for predicting the 
measured plate-mode velocities was also investigated. CPT 
was successfully used to predict the velocity of the exten- 
sional plate mode using elastic moduli calculated from lami- 
nated plate theory. The experimentally measured flexural 
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| ; J(i. 1. Tv picul Ah sigiuil in thin plate identifying extensional and flexural 
plate modes. 

mode dispersion curves were also compared with theoretical 
predictions based on CPT. The lack of agreement between 
theory and experiment at the higher frequencies demon- 
strated the limitations of C PT for composite materials. These 
are caused by the effects Of shear deformation and rotatory 
inertia which are neglected in CPT Mindlin plate theory 
which includes these effects was then used to predict the 
dispersion behavior for this mode. These predictions were in 
much better agreement w ith the experimental flexural mode 
dispersion measurements. 


These equations predict two modes of propagation. In 
general, one mode is quasi-extensional with the largest com- 
ponent of its particle displacement along the propagation di- 
rection and the other is quasi-in-plane shear w ith the largest 
component of particle displacement perpendicular to the di- 
rection of propagation in the plane of the plate. Along sym- 
metry directions, these modes are pure mode extensional and 
shear-horizontal (SHI plate modes. The dispersion relations 
for these two waves are determined by substituting a plane 
wave displacement into the equations of motion and solving 
the resulting ChristoffcTs equation. For propagation along 
the a axis or (I deg direction, this procedure produces an 
extensional mode that is pure mode and whose velocity (<*,.) 
is given by 

<Y = \A i \ iph- ^ 

Along the v axis, the extensional mode is also pure mode and 
its velocity is given by 

r ( .= \A Wp/f ■ 

For propagation at 45 deg between x and y. the wave is 
quasi-extensional with a velocity given by 


{ A 1 1 + 2 A 4- A - A + \R 


It -iA n 4 2,4 , 


4 M 1 1 4- A m* M /l :: + A (> ,, ) 


I. THEORY 

In this research, two theoretical approaches were used 
for predicting the behavior of plate mode waves in gr/ep 
composite laminates. The first theoretical predictions were 
based on CPT This is a widely used approximate theory for 
describing motion in thin plates where the wavelength (X) is 
much larger than the plate thickness (/?)• A number of au- 
thors have presented CPT in detail including Graft 11 who 
derives the equations of motion for isotropic materials and 
Whitney 12 who includes the effects of anisotropy. For exten- 
sional mode waves in a symmetric orthotropic laminate, the 
predictions of the velocity are derived from the equations of 
motion for in-plane displacements which are given by 

A 1 1 ~ + A M> 2 T ( 4 1 2 T A ) — — — ph 2 

fix ~ dx <h\ dx (ft 

(I) 

and 

d 2 v {] d~v {) d~u t) 

A ->2 + A Ml — -T- 4- { A , 2 + A 

f/y dx fix rty 

(t 2 Vn 

= ph— <2» 

p dr 

In the previous equations, n ( , and are the midplane dis- 
placements along the x and y axes which are orthogonal axes 
in the plane of the plate, p is the density, and the A if 's are the 
anisotropic in-plane stiffness coefficients obtained from 
laminated plate theory as described by Whitney 12 or Tsai and 


For flexural waves, C PT was again used to predict the 
dispersion. In this case, the plate is assumed to be under a 
stale of pure bending in which plane sections of the plate 
remain plane and perpendicular to the midplane of the plate. 
Thus, shear deformation is not included in this theory. A state 
of plane stress is assumed and the effects of rotatory inertia 
are also neglected. The CPT equation of motion for an ortho- 
tropic composite laminate in the absence of body forces is 


11 Av* 


4/7,,. + 2</J i: + 2/J„„> — T— 

<ix dx dx d\ 


4- 41), u i + t-t +/>/* -n- ^0, (7) 

dx dy ' ““ dy A ' dr 

where the D n ' s are the anisotropic bending stiffness coeffi- 
cients obtained from laminated plate theory as described by 
Whitney 12 or Tsai and Hahn. n In the previous equation, u is 
the displacement along the z axis which is normal to the 
plane of the plate. 

The dispersion behavior for the flexural mode using 
CPT is also obtained by substituting the displacement for a 
plane wave propagating along the direction of interest into 
the equation of motion. After substitution and reduction of 
terms, the resulting C PT dispersion relation along the x axis 


where c f is the velocity of the flexural mode and a> is the 
angular frequency. Along the y axis, the dispersion relation is 
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and at 45 deg the dispersion relation is 




.,'v/l 


( !) \ | + 4/) j ( , + 2( A) 1 2 + 2A) hh ) + 4/J^ + />2 1 1 


/>// 


X vw. (101 

Unlike the predictions for the extensional mode, the flexural 
mode is predicted to have a velocity that is dispersive. The 
velocity is predicted to increase as the square root of the 
frequency. 

Although CPT yielded good results for the case of the 
extensional mode, it was found to be of limited value for the 
flexural mode because of the effects of shear deformation 
and rotatory inertia. Thus, a second theory was also used to 
predict the dispersion of the flexural mode. This theory in- 
cludes the effects of shear deformation and rotatory inertia. It 
was put forth by Tang ct ai 14 following earlier work by Yang 
ct <//. ,s which was an extension of Mindlin plate theory. 10 A 
detailed presentation of this theory is beyond the scope of 
this paper. However, the dispersion behavior for a symmetric 
orthotropic laminate is obtained when the determinant of the 
following matrix of coefficients is set equal to zero 


Mu 

AY 

M 

M 2 \ 

M 21 

M 


M \ 2 

M 


where 


M\\=l>\ ] k 2 l;+2I) u ,k 2 l x l Y + l) M jL 2 ll + A"-lu> 2 . 

( 12 ) 

+ + (13) 

A/| , = /7\ 55 A7 v , < 14) 

M 2 | ~ /) | h k~ + ( /) 1 2 + )k~l J v , ( 15 ) 

(16) 


M : i-/A 44 A/ v , (17) 

= -iA 5fi kt x , UK) 

M }2 = -M 44 A/ v , 119) 


and 

M u-A + p* w 2 . (20) 

In the previous equations, l x and l v are the direction cosines 
for the propagation direction of interest, k is the wave num- 
ber. 


(P 


*./>= f 

J: hi! 


p( I 


( 21 ) 


and 


f: 1,2 

Ajj = KjKj I <(>„>/'/’■ tor /,y = 4.5. 

J : fa 1 2 


( 22 ) 


'lAhl I I. Lamina properties of AS4/3502 graphite epoxy. 


Lamina thickness 1.413 v 10 1 m 

Densilv 1550 kg/ni 1 Fiber volume 60 f r 

CL, 145,5 CiPa CL, "2.01 (JPa CL, " l L6*J (iP;i CL, 5.07 (iPa 


III Eq. (22), the k, are shear correction factors which were 
determined to yield the best agreement with three- 
dimensional elasticity theory when *7 = 5/6. The subscript / 
refers to the /th layer of the laminate and the (J,; are the 
stiffnesses for the /th layer. Solving the determinant for the 
wave number as a function of oj yields a cubic in k~. Only 
the root which approaches zero as the frequency approaches 
zero is the correct root. Once k as a function of to is known, 
the phase velocity is determined as a function of frequency 
using the relation 

Cf=-wik. (23) 


II. EXPERIMENT 

The composite laminates used in this study were made 
of AS4/3502 graphite/epoxy. All four laminates consisted of 
sixteen plys and had a nominal thickness of 2.26 mm. The 
dimensions were 0.508 m along the x direction (0 deg ply 
direction) and 0.3X1 m along the y direction. Measurements 
were made along the 0 deg tx direction), 45, and 90 deg 
directions for all four laminates. The nominal lamina prop- 
erties for this material as obtained from the manufacturer are 
given in Table I. These values were used in the laminated 
plate theory calculations to obtain the in-planu and bending 
stiffness coefficients needed for the theoretical dispersion 
calculations. 

For measurements of the extensional and flexural veloci- 
ties, a pencil lead bread was used to excite the acoustic 
waves. The waves were detected with Pannmetrics 3.5-MHz 
broadband ultrasonic transducers, These sensors have been 
shown previously by Prosser 17 to provide flat frequency, dis- 
placement sensitivity response to these low-frequency plate 
waves. The detected signals were preamplified by 40 dB with 
model I220A preamplifiers from Physical Acoustic C orpora- 
tion (PAC). There were no fillers used in the preamps. The 
amplified signals were digitized at a sampling frequency of 5 
MHz for the extensional measurements and I MHz for the 
flexural measurements with a LeCroy 6X10 transient re- 
corder, and then stored on a personal computer for analysis. 
A high sensitivity resonant AE transducer (PAC model RI5) 
was positioned next to the lead break and was used to trigger 
the transient recorder. 

For the extensional measurements, a single receiving 
sensor was used to detect the waves along the propagation 
direction of interest. The simulated AE signal was repeated 
five times while the receiving transducer was moved over a 
range of distances from the source which varied from 7.62 to 
17.78 cm in 2.54-cm increments. The first arrival time of the 
extensional mode was then determined for each of the sig- 
nals from the digitized time records and plotted as a function 
of distance from the source. A linear least-squares fit was 
used to determine the slope of the curve which is the velocity 
of the extensional mode. Phis process was repeated for each 
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of the three directions in the four plates. A typical plot of the 
extensional arrival time versus distance ot propagation is 
shown in Fig. 2. 

A Fourier phase technique was used tor the measure- 
ment of the flexural mode dispersion. This technique has 
been described by a number of authors including Sachse and 
Pao, ls Pao and Sachse, 14 Veidt and Savin' Dean. 10 and Al- 
leyne and C awley. 20 In this technique, the elastic wave is 
detected at two different distances away from the source of 
the wave along the direction of propagation of interest. The 
phase ((f) of the wave at each position at a given frequency 
(/') is determined by performing a Fourier transform on the 
signals. The phase must be unwrapped to remove the 2/?7t 
uncertainty. The phase difference (A <f) over the distance be- 
tween the two transducers (Aa) is then computed tor each 
frequency. The wave number and velocity are then calculated 
at each frequency by 


k( f) 


A tpif) 
Av 


(24) 


and 

2 77 f 2 rr / A v 
c{J) = k[f) = A<p</> ' 


(25) 


The experimental setup used for these experiments is 
shown in Fig. 3. Transducer separations of 1.91, 2.34, 3.18, 
3.81, and 4.45 cm were used. An average velocity and stan- 
dard deviation for the five different measurements was com- 
puted. The source and receivers were kept as nearly in the 



H(i. 3. Experimental setup for flexural velocity measurements in composite 
plates. 


center of the plate as possible to minimize reflections. 

Prior to computing the fast Fourier transform (FIT) to 
determine the phase, the higher-frequency extensional mode 
and the reflections arriving later in the flexural mode were 
zeroed out in the computer. Previous Fourier analysis of the 
flexural mode signals when digitized at much higher sam- 
pling frequencies (100 MHz) showed that the maximum f re - 
queney component in the flexural mode was about 200 kHz. 
Thus, aliasing was not a concern even at the low-sampling 
frequency of I MHz. 

III. RESULTS AND DISCUSSION 

The measured and theoretical extensional velocities for 
all four laminates are shown in Table II. The measured ve- 
locities were consistently lower than the predicted velocities 
with the exceptions of the 90-deg measurements in the [0 jh | 
and the (0 4 ,90 4 ] s plates and the 45-deg measurement in the 
[0J laminate. With these exceptions, the measured veloci- 
ties were all in the range from 2% to 9% less than the theo- 
retical velocities. This would seem to indicate that the manu- 
facturer s properties used in the theoretical calculations were 
somewhat higher than those in the actual material. Variations 
in material properties are common for these materials and are 
due to fiber volume variations, cure processing variations, 
and variations in resin chemistry. The 45-deg [f ) measure- 
ment was much less than theoretically predicted. It is not 
known why this discrepancy occurred and further investiga- 
tion is needed. The remaining two exceptions, propagation at 
90 deg in the [0 U J and the [0 4 ,9() 4 | s laminates, were only 
slightly higher than predicted values. 

The average measured flexural velocities for the ()-, 45-, 
and 90-deg directions in the [0,,J graphite epoxy plate are 
plotted in Figs. 4-h with the standard deviation of the mea- 
sured values indicated by error bars. The predicted velocity 
dispersion curves for CPT and Mindlin plate theory are also 
shown in these plots. The agreement between measurement 
and Mindlin plate theory is excellent for the 90-deg propa- 
gation direction. For propagation at 45 and 0 deg, the mea- 
sured values are consistently less than those predicted by 
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I 1(1. 4 Measured and theoretical llexural dispersion lor 0-dcg propagation 
in |l), (] | graph i I e/e poxy plate. 

theory. The discrepancy at 0 and 45 deg is consistent with the 
differences between theory and experiment observed for the 
extensional mode. 

The effect of shear and rotatory inertia is clear when 
CPT and Mindlin plate theory are compared in these plots. 
They are in agreement at very low frequencies in all cases 
where the approximations of ( P I are valid. The discrepancy 
increases with increasing frequency as the velocity predicted 
by ( PT increases without bound. 

It is also apparent that the difference between Mindlin 
plate theory and ( PT is much greater for the 0- and 45-deg 
directions than for the 40-deg direction. This is expected 
since the shear modulus is much smaller in comparison to 
the Young’s modulus in those directions. ( PT is based on the 
assumption of no shear deformation which implies an infinite 
shear modulus. Thus, better agreement is provided by ( PT 
w hen the ratio of the shear modulus to the Young’s modulus 
is larger 

A plot of the measured velocities and the Mindlin plate 
theory predictions for the [0,W) 4s plate is shown in Pig. 7. 
In order to better view' this complicated graph, the CPT pre- 
dictions and the experimental uncertainties are not shown in 
the plot. In this plate, the measured velocities were less than 
theoretical predictions for all three directions of propagation. 


<£ Me ri =; u r e <1 velocity 



40 80 120 160 
Frequency (kHz) 

11(1. 5. Measured and theoretical flexural dispersion tor 45-deg propagation 
in |0J graphite epoxy plate. 
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I 1(1. (>. Measured and theoretical llexural dispersion lor Wt-dcg propagation 
in |U t J graphite 'epoxy plate. 

However, the measured and theoretical velocities occurred in 
the same order with the 0-deg velocity being the largest, 
followed next by the 00-deg velocity, and with the 45-deg 
velocity the smallest. The results for the other two laminates 
w'erc similar with the measured velocities consistently less 
than predicted by theory. This again seems to indicate that 
the actual material properties are less than the nominal prop- 
erties used in the theoretical calculations and is consistent 
with the extensional measurements. 

In summary, t he velocity of the extensional plate mode 
and the dispersion of the llexural plate mode was measured 
along three directions of propagation in four graphite epoxv 
laminates. Theoretical predictions of the dispersion of these 
inodes were made with CPT and Mindlin plate theory. CPT 
yielded good agreement with the extensional velocity but 
was shown to be of limited value in predicting the dispersion 
of the flexural mode because it assumes that the effects of 
shear deformation and rotatory inertia are negligible. Mind- 
lin plate theory, w hich includes these effects, gave much bet- 
ter agreement with the measured flexural dispersion. How- 
ever, there was a consistent discrepancy between theory and 
experiment believed to be due to variations in actual material 
properties from those used in the calculations. 


1500 ~ 



I 1(1. 7. Measured and theoretical flexural dispersion lor U-. 45-. and 90-dcg 
propagation directions in |U, t )()| 4% graphite/epoxy plate. 
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